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Abstract 



We introduce a class of exactly solvable boson models. We give explicit 
analytic expressions for energy eigenvalues and eigenvectors for an sd-boson 
Hamiltonian, which is related to the SO (6) chain of the Interacting Boson 
Model Hamiltonian. 



I. INTRODUCTION 

A number of complex quantum systems can be described algebraically using bosonic 
degrees of freedom. For example even-even nuclei can be described in terms of a system of 
one boson with angular momentum L = 0, called s, and five components of another boson 
with angular momentum L = 2, called d^, fi = 0,±1,±2. This model, called Interacting 
Boson Model, has been very successful in describing the low-lying collective states of even- 
even nuclei [1,2] A similar boson model, using an angular momentum L = 1 (p) boson in 
addition to the s-boson has been used to describe the structure of molecules [3] and hadrons 
[4]. The s-d boson model has an SU(6) symmetry and the s-p boson model has an 577(4) 
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symmetry. Analytic solutions of these models are given whenever the Hamiltonians can be 
written in terms of the Casimir operators of the SU(6) (or SU(4)) group and its subgroups. 
In particular for the 5*0(6) chain the Hamiltonian contains the boson pairing operator. 

Since the boson creation operators are spherical tensors and the annihilation operators 
can be written as spherical tensors after the introduction of a suitable phase (i.e. = 
(— lyd-n) we will use the dot product notation of spherical tensors, (•), below. It is easy to 
show that the three operators 

S + = ±(<t-<ft), (L1) 



S- = \(d-S), (1.2) 

and 

S ° = \j2(dtd, + d,dl), (1.3) 

generate the d-boson quasi-spin group SU(1, 1). This group can also be represented using 
only the s- or p-bosons. The eigenvalues of the operator 5° are given by (n/2) + (1 + a)/4, 
where n is the total number of the bosons and the quantity a is 0, 1, and 2 for the s-,p-, 
and d-bosons, respectively. For the 50(6) symmetry chain one needs to include both the s- 
and d- bosons in the realization of the algebra: 

5+ = I[(dt.rft)-( a t a t )]> (1 . 4) 



S- = -[(d ■ d) - (ss)], 



;i. 5 ) 



and 



L V 



(1.6) 



The pairing operator of the 50(6) limit is then P 6 = 5 + 5 , where 5 ± are given by Eqs. 
(1.4) and (1.5). 

Our goal in this paper is to consider a system of s- and d-bosons near, but not necessarily 
at the 50(6) limit. It is worth to point out that similar efforts were undertaken by others. In 
Ref. [5] the U(5) to SU (3) transition in the interacting boson model was considered. In Ref. 
[6] an exactly solvable model of bosons with a repulsive pairing and single-particle energies 
proportional to the angular momentum quantum number was introduced. This model was 
further elaborated in Ref. [7]. In these models the eigenvalues and the eigenfunctions are 
written down in terms of variables that need to be numerically evaluated by solving non- 
linear equations. Our aim here is to introduce a model for which an explicit analytic energy 
level expression can be given for all levels and analytic expressions can be given for at least 
few lowest eigenstates. A preliminary version of this work was presented in Ref. [8]. 
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II. THE MODEL 



To present our model we first introduce the operator 



* + (A) = ~ 2 



1-A 



A 



(2.1) 



Clearly S + (\) is proportional to the operator in Eq. (1.4) when the parameter A = 1/2. 
However it represents a general mixing of s and d bosons when A ^ 1/2. Similarly we 
introduce the operators 



s-(A) = 5 



(d ■ d) ss 
1-A ~ X 



and 



S°(A) 



^luid^dn + d^djj) 



(s^s + ss^) 



1-A A 

In this paper we show that the eigenvalues and the eigenstates of the Hamiltonian 



(2.2) 



(2.3) 



H = 



1 



N d 



A A 



+ 



N d Ns 
- A) 2 A 2 



-W{\) 



N d 



N, 



A A 



+ S + (\)S-(X), (2.4) 



can be analytically calculated. Here N d and N s are the d- and s-boson number operators, 
and the function W(X) is defined below. 

The operators given in Eqs. (1.1), (1.2), and (1.3) and their s-boson counterparts satisfy 
the SU(1, 1) commutation relations 



[S^~, Sj 



2SijSj, 



and 



[Si, Sf] — ±5ijSf, 



(2.5) 



(2.6) 



where etc. represent different SU(1, 1) algebras written using only the s,d, etc. bosons. 
Using Eqs. (2.5) and (2.6) we obtain the commutation relations 



g°(A)-gy) 

A — ji 



(2.7) 



[s%x),sH^] = ± S±{X l S±i(i \ 

A — /! 



(2.8) 



[5 (A),5°(^)] = [5 ± (A),5 ± (|i)] = 0. 



(2.9) 



Eqs. (2.8) were first studied by Gaudin [9] in the context of interacting spins on a lattice. 
The limit of these commutators when X — /i can easily shown to be 
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[S+(\),S-(\)] = -2 5 -^, (2.10) 

[S°(A),^(A)] = ±^^, (2.11) 
Next it is easy to show that the parametric Hamiltonians 

H(X) = S°(X)S°(X) - V(A)^T(A) - V(A)S+(A) (2.12) 
form a one-parameter family of mutually commuting operators: 

[H(\),H(jm)] = 0. (2.13) 

Starting from a lowest weight vector and using J + (X) as step operators, one can diagonalize 
them simultaneously. Lowest weight vector |0 > is the boson vacuum: 

S-(A)|0>=0. (2.14) 

Note that our choice of the ground state is different from what is normally taken in the 
SO(6) limit, which has the maximum population of the s-bosons. The function W(X) we 
introduced earlier in the Hamiltonian of Eq. (2.4) is obtained by the action of the operator 
S°(\) on the vacuum state: 

S°{X)\0>=W{\)\0> . (2.15) 

For the s-d boson model W(X) takes the form 

1 



W(X) { 



(2.16) 



.1-A A. 

The boson vacuum itself is an eigenvector of H(X): 

H(X)\0>=e (X)\0> (2.17) 

with the eigenvalue 

e (A) = W{X) 2 - W\\). (2.18) 
It can be shown that a vector of the form 

|£>=S+(0|O> (2.19) 
is also an eigenvector of H(X) if £ the solution of 

W(0 = 0. (2.20) 

In this case the eigenvalue is given by 
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Cl (A) = 6 (A) - 2^. (2.21) 

Using the algebra given in Eqs. (2.8) it is easy to show that the state 

|6, >= S+i&S+fa) . . . S + (£ n ) |0 > (2.22) 

is an eigenvector of H(X) if the quantities £1, £ 2 , • • ■ satisfy the following system of equations: 

n 1 

W((«) = Er77 for « = l,2,...,n. (2.23) 
The eigenvalue of the Hamiltonian in Eq. (2.12) on the state of Eq. (2.22) is given by 

^H, W - 2 e!". (2.24) 

a=l ^ Q 

Comparing the Hamiltonians of Eq. (2.4) and (2.12) and using Eqs. (2.11), it follows that 
the energy eigenvalues of the Hamiltonian TC are given by 

kw = 2± m^M . (2 . 25) 

All the above considerations are equally applicable to both real and complex values of A. 
However only when A is real then our Hamiltonian is Hermitian with real eigenvalues. This 
is true even when the solutions of Eq. (2.23) are complex. 



III. CALCULATION OF THE EIGENVALUES AND THE EIGENFUNCTIONS 

In order to evaluate energy eigenvalues and eigenfunctions we will use the general ex- 
pression 



W(X) 



A 



B 



a\ — A a 2 — A 



(3.1) 



so that our calculations are equally applicable to the mixture of other bosons. Inserting Eq. 
(3.1) into Eq. (2.25) it is easy to show that 



E n (\) = 2 £ 



+ 



B 



(ai - A) (ai - £ Q ) (a 2 - A) (a 2 - C 



(3.2) 



It turns out that an expression for energy eigenvalues can be obtained without explicitly 
calculating the solutions of Eq. (2.23). From Eq. (2.23) it is fairly straightforward to show 
that 



E W(£a) = 0. 

Q=l 



(3.3) 



Similarly considering the expression 
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E taW(t a ) = £ 



a=l 



a=l 



E 



Co 



+ E 



Co 



0< 



a Ca £/3 /3 >a Ca £/3 



and interchanging the dummy indices a and (3 in the last sum one finds that 

a=l (a,/3)pairs ^ a 



Using Eqs. (3.1) and (3.3) one then gets 

n 

A E : 

Inserting Eq. (3.1) into Eq. (3.5) and employing the identity 



- 1 



we get 



^E 



i +^E 



a 2 



n(n — 1) 



Using Eq. (3.6) in Eq. (3.8) one then obtains 



B(a 2 - a x ) E 



n(n — 1) 



+ (A + B)n. 



Using Eq. (3.6), the energy eigenvalues of Eq. (3.2) take the form 



E n (X) = 2B 



E 



,a 2 - A ai - A/ ^ a 2 - £ a 
Finally inserting Eq. (3.9) into Eq. (3.10) we obtain 



E n (X) 



-2 



n(n — 1) 



+ (A + B)n 



(ai - A)(a 2 - A) 

In our model we took a x = 1, a 2 = 0, A = 5/4, and 5 = 1/4 resulting in 

1 



E n (\) 



A(l-A) 



[n 2 + 2n]. 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



Even though we did not need explicit solutions of Eq. (2.23) to evaluate the energy 
eigenstates, we nevertheless need to find these solutions to be able to write down the energy 
eigenstates. To achive this goal we first introduce a change of variables 



ia = a 2 + ( a (ai - a 2 ). 
Using these new variables the Bethe ansatz equation Eq. (2.23) takes the form 



(3.13) 
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E 



B 



A 



a Ca C/3 Ca 1 Ca 



for a — 1,2, ... ,n, 



(3.14) 



where B = 1/4 and A = 5/4 in the s-d boson problem. In Ref. [10], Stieltjes had shown 
that the polynomial 



Pr,(0=II(C-Ca) 
a=l 

satisfies the hypergeometric differential equation 

x(l - x)p£(x) + [2B - 2(A + B)]p' n (x) + n(n + 2A + 2B- l)p n {x) = 0. 
Hence ( a are the roots of the polynomials 



E 

k 



(-n) k (n + 2A + 2B-l) k k 
(2B) k k\ 



where 



(a) n = a(a + 1) • • • (a + n — 1) 
for n — 1, 2, • • • and (a) = 0. So for example the first excited eigenstate is 



A + B 



|0>, 



and the second eigenstate is 

(dt • rft) 



(A + 5 + l)^ 



S t s t 



+ 1 - 7 25 + 1+7 



S t s t 



2^1 + 1 + 7 ~ 25+1-7 



where 
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, (2A+1)(25 + 1) 
2A + 2B + 1 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



|0 >, (3.20) 



(3.21) 



Finding the eigenstates of the third, fourth, etc. excited states simply requires finding the 
roots of a cubic, quartic, etc. polynomial. Even though for the third and fourth order 
polynomials this can be done analytically, the resulting expressions are not particularly 
illuminating. Of course the roots can be calculated numerically for any order. 
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